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ABSTRACT 

Previous observations of the middle-aged pulsar Geminga with XMM-Newton and Chandra have 
shown an unusual pulsar wind nebula (PWN), with a 20" long central (axial) tail directed opposite 
to the pulsar's proper motion and two 2' long, bent lateral (outer) tails. Here we report on a deeper 
Chandra observation (78 ks exposure) and a few additional XMM-Newton observations of the Geminga 
PWN. The new Chandra observation has shown that the axial tail, which includes up to three brighter 
blobs, extends at least 50" (i.e., 0. 06^250 pc) from the pulsar (CZ250 is the distance scaled to 250 pc). 
It also allowed us to image the patchy outer tails and the emission in the immediate vicinity of 
the pulsar with high resolution. The PWN luminosity, £o.3-8k c v ~ 3 x 10 29 g?25o er § s_1 > ^ s lower 
than the pulsar's magnetospheric luminosity by a factor of 10. The spectra of the PWN elements 
are rather hard (photon index r ~ 1). Comparing the two Chandra images, we found evidence of 
PWN variability, including possible motion of the blobs along the axial tail. The X-ray PWN is the 
synchrotron radiation from relativistic particles of the pulsar wind; its morphology is connected with 
the supersonic motion of Geminga. We speculate that the outer tails are either (1) a sky projection of 
the limb-brightened boundary of a shell formed in the region of contact discontinuity, where the wind 
bulk flow is decelerated by shear instability, or (2) polar outflows from the pulsar bent by the ram 
pressure from the ISM. In the former case, the axial tail may be a jet emanating along the pulsar's 
spin axis, perhaps aligned with the direction of motion. In the latter case, the axial tail may be the 
shocked pulsar wind collimated by the ram pressure. 

Subject headings: pulsars: individual (Geminga) — stars: neutron — stars: winds, outflows — X-rays: 
stars 



1. INTRODUCTION 

Pulsars lose their spin energy via relativistic pulsar 
winds (PWs) of charged particles. The PW shocks in 
the ambient medium and forms a pulsar wind nebula 
(PWN) whose synchrotron radiation can be observed in 
a very broad energy range, from the radio to TeV 7- 
rays (see Kaspi et al. 2006, Gaensler & Slane 2006, and 
Kargaltsev & Pavlov 2008 [KP08 hereafter] for recent re- 
views). The shocked PW is confined between the termi- 
nation shock (TS) and contact discontinuity (CD) sur- 
face that separates the shocked PW from the shocked 
ambient medium between the CD and the forward shock 
(FS). The shapes of the TS, CD, and FS depend on 
the wind outflow geometry and the ratio of the pulsar's 
speed to the sound speed in the ambient medium (the 
Mach number), A4 — v psr /c s . In particular, if the pul- 
sar moves with a supersonic speed, M. 3> 1, and the 
preshock PW is isotropic, then the TS, CD, and FS ac- 
quire bow-like shapes ahead of the pulsar, with the TS 
apex ("head") at a distance Pts.Ii ~ (E/4Trcp Tairi ) 1 ^ 2 , 
where p lam = Pamb^psr i s the ram pressure, p am b the den- 
sity of the ambient medium (e.g., Bucciantini et al. 2005; 
hereafter B+05 ). The shocked PW forms a tail behind 
the pulsar, with a flow speed significantly exceeding the 
pulsar's speed (Romanova et al. 2005; B+05). 

Among ~60 PWNe detected by Chandra, about 20 
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PWNe show such bowshock-tail morphologies (KP08). 
Such tails have been observed, for instance, behind the 
pulsars J1747-2958 (Gaensler et al. 2004), J1509-5850 
(Kargaltsev et al. 2008), B0355+54 (McGowan et al. 
2006), and B1929+10 (Misanovic et al. 2008), with very 
different spindown ages, r sd = P/(2P) = 26, 160, 620, 
and 3100 kyr, respectively. We should note, however, 
that the detailed shape of the detected bowshock-tail 
PWNe is often different from the idealized models, es- 
pecially in the immediate vicinity of the pulsar, possibly 
because of anisotropics of the pulsar outflows. For in- 
stance, by analogy with a few bright, well-resolved PWNe 
around young pulsars moving with subsonic velocities 
(such as the Crab PWN; Weisskopf et al. 2000), one 
can expect that the pulsar outflows consist of equato- 
rial and axial components, with respect to the spin axis, 
which are responsible for the "tori" and "jets" observed 
in these torus-jet PWNe (KP08). 

One of the most peculiar PWNe has been detected 
around the famous Geminga pulsar (PSR J0633+1746). 
Geminga was discovered as a 7-ray source 7195+5, with 
the SAS-2 satellite (e.g., Thompson et al. 1977). The 
period of Geminga, P = 237 ms, was discovered by 
Halpern & Holt (1992) in X-ray observations with the 
Rontgen Satellit (ROSAT), and the period derivative, 
P = 1.1 x 10~ 14 s s _1 , was first measured by Bertsch 
et al. (1992) in 7-rays with the Compton Gamma Ray 
Observatory (CGRO). The period and its derivative cor- 
respond to the spindown age t s( j = 340 kyr and spin- 
down power E — 3.3 x 10 34 erg s _1 . The Geminga 
pulsar has also been detected in the optical (Halpern 
& Tytler 1988; Bignami et al. 1988), near-IR (Kopt- 
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sevich et al. 2001), and UV (Kargaltsev et al. 2005). 
The distance to Geminga, d = 250^^° pc, was esti- 
mated from its annual parallax measured in observa- 
tions with the Hubble Space Telescope (Faherty et al. 
2007). Its proper motion, 178.2 ± 1.8 mas/yr, corre- 
sponds to the transverse velocity, vj_ ~ 211(i25o km s _1 
[where g?250 = d/(250pc)]. As this velocity consider- 
ably exceeds the typical sound speed in the interstellar 
medium (ISM), c s ~ 10-30 km s _1 , one should expect 
Geminga to be accompanied by a bowshock-tail PWN, 
with i?TS,h = 1.1 x 10 16 n _1 / 2 (d25o/ sini)^ 1 cm, which 
corresponds to ~ 2.9"n~ 1 / 2 (d 250 / sini) -2 , where i is the 
angle between the pulsar's velocity and the line of sight, 
and n = p^b/m^. 

XMM-Newton observations of Geminga in 2002 April, 
reported by Caraveo et al. (2003; hereafter C+03), re- 
vealed two w 2' long tails behind the pulsar, approxi- 
mately symmetric with respect to the sky projection of 
the pulsar's trajectory (see Fig. 1), with a luminosity of 
- 10 29 erg s" 1 in the 0.3-5 keV band. C+03 suggested 
that these tails are associated with a bowshock generated 
by the pulsar's motion, and, using the one-zone bowshock 
model by Wilkin (1996) 4 , predicted that the head of the 
bowshock, 20"-30" ahead of the pulsar, is hidden in the 
bright wings of the pulsar point spread function (PSF) 
in the XMM-Newton image. 

The Geminga field was observed in 2004 (Sanwal et 
al. 2004; Pavlov et al. 2006 [hereafter P+06]) with the 
Chandra Advanced CCD Imaging Spectrometer (ACIS), 
whose resolution, « Of! 5, is much better than that of the 
XMM-Newton detectors. The most interesting finding 
of that observation was the detection of an axial tail be- 
hind the pulsar aligned with the direction of the pulsar's 
proper motion (P+06; de Luca et al. 2006; see Fig. 2, 
top). The axial tail, with a luminosity ~ 2 x lO 29 ^ 2 ^ 
erg s _1 , was seen up to 25" from the pulsar, almost up 
to the boundary of the held of view (FOV). P+06 sug- 
gested that the axial tail may be a jet emanating from 
the pulsar magnetosphere. In addition to the axial tail, 
a faint arc-like structure was detected 5" -7" ahead of 
the pulsar (but no emission at 20"-30", contrary to the 
C+03 prediction), and a 3 a enhancement, apparently 
connecting the arc with one of the outer tails (south of 
the axial tail), was noticed (P+06). No emission was de- 
tected from the other (northern) outer tail in that short, 
« 20 ks, exposure. 

To image the whole extent of the Geminga PWN and 
study its tails in more detail, we observed this field with 
Chandra ACIS in 2007, with a longer exposure and a 
larger FOV. In this paper, we report the results of this 
observation and compare them with the previous find- 
ings. We describe the data analysis and the observational 
results in § 2, and discuss the implications of these results 
in §3. 



4 The one-zone model assumes an instant mixing of the PW 
matter with the ambient matter at the shock, so that there is no 
CD, and the TS and FS suraces coincide with each other. The 
numerical bowshock models (e.g., Bucciantini 2002) have shown 
that the shape of the one-zone shock is approximately similar to 
that of the FS but very different from the TS shape. 



2. DATA ANALYSIS AND RESULTS 
2.1. Observations and data reduction 

The Geminga field was observed with Chandra ACIS 
on 2007 August 27 for 78.12 ks (ObsID 7592). The 
observation was taken in Timed Exposure (TE) mode, 
with the frame time of 3.24 s. After removing 20 s of 
high background and correcting for the detector dead 
time, the scientific exposure time (live time) is 77,077 
s. To maximize the signal-to-noise ratio (S/N) for the 
very faint PWN emission, we imaged the field onto the 
front-illuminated 13 chip, which has a lower background 
than the commonly used (and slightly more sensitive) 
back-illuminated S3 chip 5 . We used the Very Faint 
telemetry format to provide a better screening of back- 
ground events 6 . As putting the target at the ACIS-I 
aimpoint (near the corner of the 13 chip) could result 
in chip gaps crossing the PWN image, we moved the 
focus to the middle of node 2 on the 13 chip (SIM-Z 
= -7.42 mm = -2.'53) and applied the AY = -1.'6 off- 
set to put the pulsar at > 2' from the chip boundaries. 

To obtain a deeper PWN image and examine a possi- 
ble PWN variability, we also used the previous Chandra 
observation of Geminga carried out on 2004 February 7 
(ObsID 4674; 18,793 s scientific exposure ). The details 
of that observation have been described by P+06. Here 
we only mention that the observation was taken in Faint 
telemetry format, and the Geminga pulsar and its PWN 
were imaged on 1/8 subarray of the S3 chip (w V x 8' 
FOV), which reduced the pileup in the pulsar image but 
did not allow us to image the whole PWN. 

We have used the Chandra Interactive Analysis of Ob- 
servations (CIAO) software (ver. 4.0; CALDB ver. 3.4.0) 
for the ACIS data analysis, starting from the level 1 event 
files. We have applied the standard grade filtering and 
used the energy range 0.3-8 keV to minimize the back- 
ground contribution. We have also applied the exposure 
map correction, but found that the effects of nonuniform 
exposure and nonuniform CCD response in the PWN re- 
gion are small (except for the boundary of the 1/8 sub- 
array used in the observation of 2004). 

To confront the high-resolution Chandra data with 
the XMM-Newton results, we also used the data ob- 
tained with the MOS1 and MOS2 detectors of the Euro- 
pean Photon Imaging Camera (EPIC) on board XMM- 
Newton. In addition to the observation of 2002 April 
4-5 (ObsID 011117010; 77.97 ks scientific exposure, af- 
ter removing the periods of high background) reported by 
C+03 and P+06, we also used the data sets obtained in 
observations 0201350101 of 2004 March 13 (16.23 ks sci- 
entific exposure), 031159100 of 2006 March 17 (4.46 ks), 
0400260201, 2006 October 2 (19.70 ks), and 0400260301 
of 2007 March 11 (23.83 ks). The total effective expo- 
sure of the five observations is 142.18 ks. All the ob- 
servations were taken with medium filter in Full Frame 
mode, providing a 30' diameter FOV. The data reduction 
was performed with the Scientific Analysis System (SAS) 
package (ver. 8.0.0). Good events with patterns 0-12 and 
energies within the 0.3-8 keV range were selected for the 
data analysis. 

5 See §6.16 of The Chandra Proposers' Observatory Guide 
(POG), http:/ /asc. harvard.edu/proposcr/POG. 

6 See § 6.14 of the Chandra POG. 



New X-ray Observations of the Geminga PWN 



3 



2.2. Images and spectra of the PWN elements 

The Chandra ACIS data of 2007 provide the high- 
resolution image of the entire Geminga PWN for the first 
time (see the top panel of Fig. 2 and Fig. 3). Below we 
will describe the observed properties of the PWN ele- 
ments, and compare them with the results of the 2004 
ACIS observation. 

To calculate the net source counts N$ in the area As 
of a PWN element, we use the formula N$ — Nt — 
(As/Ab)Nb, where Nt is the total number of counts 
detected from the area As, and Nb is the number of 
background counts detected from the area As- Then 
the 1 (7 error of Ns and the signal-to-noise ratio are 
given by SN S = [N s + (1 + A s / A B ){A S / A b )Nb\ 1/2 and 
S/N = N S /SN S , respectively. 

For the analysis of the 2004 data, we use a source- 
free rectangular background region (As — 3439 arcsec 2 , 
Nb = 129 counts in the 0.3-8 keV band) to the north 
of the Geminga pulsar. For the 2007 data, we use the 
background measured from a source-free rectangular re- 
gion, with the area As — 5399 arcsec 2 , in the northeast 
portion of the ACIS-I3 chip (unless stated otherwise). 
This region contains Nb — 416 counts in the 0.3-8 keV 
band, which corresponds to the background brightness of 
1.0 x 10 -6 counts arcsec" 2 s _1 , a factor of 2 lower than 
in the Chandra observation of 2004. The values of N$ 
and S/N for the PWN elements are given in Table 1. 

For the spectral analysis, we have used the XSPEC 
package (ver. 12.4.0) and fit the spectra with the ab- 
sorbed power-law (PL) model (wabs*powerlaw), with 
the fixed hydrogen column density N-r = 1.1 x 10 20 cm~ 2 
(Halpern & Wang 1997; de Luca et al. 2005). As the 
number of counts in the PWN is small, we use the maxi- 
mum likelihood method (C-statistic) for spectral fitting. 
Table 1 provides the values of the photon index T and 
the flux F of the PWN elements. 

2.2.1. Axial tail 

The brightest feature of the Geminga PWN in the 
Chandra data of 2004 is the axial tail (A-tail hereafter) , 
seen up to at least 25" from the pulsar in the direction 
opposite to the pulsar's proper motion (see P+06 and 
Fig. 2, middle). In the image from the 2007 observation, 
we see the A-tail up to at least 50" (0.06 c^so pc) from 
the pulsar (Fig. 2, top), with 83±12 source counts within 
the region of 706 arcsec 2 area shown by the solid lines 
in Figure 3. The PL fit of its spectrum (see Table 1 and 
Fig. 4) gives the photon index T = 1.8±0.3 and the 0.3-8 
keV luminosity L — (0.9 ± 0.1) x 10 29 erg s _1 (assuming 
an isotropic emission at d = 250 pc), versus T = 1.3 ±0.3 
and L = (1.6 ± 0.3) x 10 29 erg s" 1 in the 2004 data, as 
measured in the 118 arcsec 2 area rectangle that contains 
46 ± 7 counts (shown in the middle panel of Fig. 2). 

The A-tail looks patchy in both the 2004 and 2007 ob- 
servations, with some "blobs" standing out. The blobs, 
labeled A, B, and C in Figure 2 are at the distances 
of about 9", 20", and 43" from the pulsar, respectively 
(blob A is seen in the 2004 image, while blob B and blob 
C are seen in the 2007 image). They contain 24.9 ± 5.1, 
12.8 ± 3.9, and 25.8 ± 5.3 source counts, respectively, 
within the 3"-radius circles around their centers. The 
analysis of the brightness distribution along the A-tail 
shows that the blobs are significant at > 3er levels (i.e., 



they are not just statistical fluctuations of the bright- 
ness distributions). For instance, the number of counts 
in the 2" radius circle around the center of blob B (15 
counts) exceeds the average number of counts per the 
same 12.6 arcsec 2 area in the A-tail (2.44 ± 0.21 counts) 
at the 3.2 a level. The nonuniform surface brightness 
distribution along the A-tail, and the difference of these 
distributions in the 2004 and 2007 images are shown in 
Figure 5. 

The brightest in the 2007 data is blob C at the ap- 
parent end of the tail, centered at a = 06 h 33 m 51?71, 
6 = +17°45'53'/1 (J2000). Because of the small number 
of counts, we cannot firmly determine whether the blob 
corresponds to a point source or an extended one. In- 
terestingly, the end portion of the tail looks attached to 
this blob, while the tail looks detached from the pulsar 
in both the 2004 and 2007 images. Therefore, one could 
even speculate that the tail might belong not to Geminga 
but to some unrelated field object (e.g., it might be a jet 
of an active galactic nucleus [AGN] , accidentally oriented 
toward Geminga in the sky projection). 

To check whether the blobs are indeed associated with 
the tail or they may be background sources, we examined 
the optical/NIR catalogs. We found no optical counter- 
parts to blob A and blob B, but we found an object at 
a = 06 h 33 m 51!69, S = +17°45'54'/2 (J2000), about 1" 
from the center of blob C (the coordinates are from the 
USNO-B1.0 catalog [Monet et al. 2003], with the quoted 
mean uncertainties of 0'.'075 and 0'.'094 in a and 6, re- 
spectively). Based on the magnitudes and colors (e.g., 
V = 17.7 ± 0.3 [GSC2.3 catalog; Lasker et al. 2008], 
J = 16.18 ± 0.09, H = 15.54 ± 0.11, K = 15.6 ± 0.2 
[2MASS catalog; Cutri et al. 2003]), this object could be 
a background K star. Such a star could contribute to 
the X-ray emission of blob C. The observed X-ray flux in 
the 3" radius aperture is -Fo.3-3.5koV ~ 3.8 x 10~ 15 erg 
cm~ 2 s" 1 . According to Maccacaro et al. (1988), the X- 
ray/optical flux ratio, log(Fo.3_3.5k e v/-PV) ~ —2.0, cor- 
responds to a K or M star, and it excludes an AGN [for 
which —1.2 < log(F .3_3.5kcv/-FV) < +1-6) as the source 
of the X-ray and optical emission (hence the tail is not 
an AGN jet). Thus, we cannot rule out the possibility 
that a K star, accidentally projected onto the A-tail, is 
at least partly responsible for the brightened end of the 
tail in the 2007 image. However, if we exclude blob C, 
the tail's luminosity decreases by a factor of 1.5, but the 
spectral slope remains virtually the same (see Table 1 
and Fig. 4). This suggests that the star's contribution 
does not dominate in the blob C emission, but, because 
of the small number of counts and large statistical errors, 
we cannot firmly conclude on the nature of blob C. 

In the observation of 2004, blob C was imaged onto an 
underexposed part of the FOV (because of the dither), 
between the dashed lines in the middle panel of Figure 2. 
Taking into account the shorter effective exposure of that 
observation (but the higher sensitivity of the S3 chip) , we 
expect 4.2±2.1 counts to be detected in the 2" radius cir- 
cle around the position of the blob C centroid; however, 
there are no counts within that circle. This may suggest 
some variability of the source, but the statistical signifi- 
cance of this difference is marginal (e.g., the probability 
of detecting zero counts when 4.2 counts are expected is 
0.0145, which corresponds to a 2.4 a significance). 

Based on the ACIS count rate of blob C in 2007, one 
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could expect to detect about 90 counts in the 15" ra- 
dius aperture in the 142 ks MOS1+MOS2 exposure, but 
we found only 49 ± 12 counts in the MOS data. Fur- 
thermore, we note that the blob C position is projected 
onto the wings of the pulsar PSF in the XMM-Newton 
images (see Fig. 1), whose contribution to the number 
of extracted counts is difficult to evaluate because of the 
"spiky" shape of the PSF. Anyway, the number of counts 
expected for blob C in the XMM-Newton data signifi- 
cantly exceeds the measured one, suggesting variability 
of blob C. 

The A-tail images (Fig. 2) look appreciably different in 
the 2004 and 2007 data, in both the overall flux and the 
surface brightness distribution (see Fig. 5). For instance, 
the flux, 2.2 x 10~ 14 erg cm -2 s -1 , in the 118 arcsec 2 
area of the tail in the 2004 data is a factor of 6 higher 
than the flux from the same area in the 2007 data (the 
difference between the count rates is significant at the 5 a 
level, with account for a factor of 1.6 higher sensitivity of 
the S3 chip compared to that of the 13 chip, for T = 1.8 
and A H = 1-1 x 10 20 cm" 2 ). 

The different positions of the blobs in the 2004 and 
2007 images suggest that the blobs are moving along the 
A-tail (perhaps similar to the blobs in the Vela pulsar 
jet; Pavlov et al. 2003). One might even speculate that, 
for instance, blob B is, in fact, blob A that had moved 
11" (4 x 10 16 <i25o cm) in 3.5 yr between the observa- 
tions (which would correspond to the transverse velocity 
of «biob,_L ~ 3700 c?25o km s -1 ). However, as blob B could 
also form independently after the disappearance of blob 
A, this will remain a speculation until the characteris- 
tic blob lifetime is estimated in a series of monitoring 
observations. 

In the large-scale Chandra image of 2007 (Fig. 3), one 
can see a possible faint extension of the A-tail (within 
the dashed polygon in Fig. 3). This faint portion has a 
factor of 1.5 higher observed flux (but a factor of 2.4 lower 
average brightness) than the bright portion (see Table 1). 
However, its statistical significance is only 1.7 c, and it 
is not seen in the deep XMM-Newton image (see Fig. 1). 
Therefore, this "faint portion" most likely represents a 
string of background fluctuations accidentally aligned in 
the A-tail direction. 

2.2.2. Outer tails 

The two "outer tails" of the Geminga PWN, seen up 
to ~ 2' from the pulsar, were originally detected in 
the XMM-Newton observations of 2002 April (C+03). 
Adding four shorter observations of 2004-2007, which 
increases the total exposure by a factor of 1.8, shows 
qualitatively the same picture (Fig. 1). The two tails 
are approximately symmetric with respect to the pul- 
sar's trajectory in the sky, forming a horseshoe-like struc- 
ture. The southern and northern tails (we will call them 
the S-tail and N-tail, for brevity) are seen up to 3.'1 
and 2.'7 from the pulsar, respectively, in the summed 
MOS1+MOS2 image. Their typical width, - 20"-30", 
is comparable to the XMM-Newton angular resolution. 
The tails in the pulsar vicinity (within « 30") are im- 
mersed in the bright pulsar's image. The spectrum of 
the combined emission from the two tails, extracted from 
the 120" x 40" elliptical regions shown in Figure 1 (about 
560 source counts), can be described by a PL model with 
T = 1.7±0.2 (xl = 1.36 for 50 degrees of freedom [d.o.f], 



for fixed Ah = 1.1 x 10 20 cm -2 ). The unabsorbed flux 
and luminosity of the two tails, are F unabs ~ 1.7 x 10~ 14 
erg cm" 2 s" 1 and L ~ 1.3 x lO 29 ^ 2 ^ erg s _1 , in the 0.3- 
8 keV band. Both the total flux and the spectral slope 
are consistent with those obtained by C+03. The con- 
tribution of the S-tail into the total energy flux is about 
76%. The average specific intensities in the S-tail and N- 
tail elliptical regions arc 3.4 x 10~ 18 and 1.1 x 10~ 18 erg 
cm -2 s -1 arcsec -2 , respectively. Fitting the spectra of 
the S-tail and N-tail separately, we obtained T = 1.4±0.3 
and 2.7 ± 1.1, respectively. The apparently large differ- 
ence between the spectral slopes, Ar = 1.3 ± 1.1, is not 
statistically significant. The spectra can also be fitted 
by the models for emission of an optically thin thermal 
plasma (e.g, kT = 8±J° keV for the fit of the S-tail + 
N-tail spectrum with the mekal model: \t = 1-36 for 50 
d.o.f.). 

Thanks to its high angular resolution, Chandra obser- 
vations make it possible to image the tails in the pulsar 
vicinity and resolve the tail structure. In the short ob- 
servation of 2004 the N-tail was not detected, while the 
initial portion of the S-tail was detected with about 3 a 
significance (P+06). The entire extent of the tails could 
not be seen because the 1/8 subarray was used. 

In the Chandra data of 2007 we have detected both 
outer tails. For the analysis, we divide each of the tails 
into two parts. The bright initial parts (up to 2' and 40" 
from the pulsar, for the S-tail and N-tail, respectively) 
are delineated by solid polygons in Figure 3, while the 
longer faint parts (up to 4.'2 and 3'6 from the pulsar, 
respectively) are shown by dashed polygons. (We note 
that the bright portion of the N-tail and a substantial 
part of the bright portion of the S-tail are hidden behind 
the pulsar image in the XMM-Newton data.) The sta- 
tistical significance of the faint parts is marginal in the 
Chandra data (see Table 1), but their reality is supported 
by the XMM-Newton data (see Fig. 6, where the XMM- 
Newton brightness contours are overlaid on the Chandra 
image) . The end parts of the S-tail in the Chandra and 
XMM-Newton images are slightly shifted with respect to 
each other, which might suggest possible variability of 
the outer tails. The Chandra image resolves the tails in 
the transverse dimension, showing > 20" widths, but the 
image is not deep enough to infer the brightness distribu- 
tion across the tails. Also, the tails do not show sharp-cut 
outer boundaries, perhaps because of the same reason. 

Figure 6 also shows the locations of four optical- NIR 
sources that are projected onto the tails. We have al- 
ready discussed one of them (#1 in Fig. 6, whose posi- 
tion coincides with blob C in the A-tail). The other three 
sources (numbered 2, 3, and 4 in Fig. 6, with V magni- 
tudes of 17.6, 15.3, and 13.4, respectively, from the NO- 
MAD catalog; Zacharias et al. 2004) are projected onto 
the S-tail. Source 3 and source 4 are likely an F star 
and a K star, respectively, based on their optical-NIR 
colors, while the magnitude errors for the fainter source 
2 are too large to determine its nature. The 2.55" ra- 
dius circles around the positions of the sources 2, 3, and 
4 contain 7, 1, and 11 counts, respectively. Even if the 
X-ray emission at these locations is due to the optical 
sources, the total number of detected source counts from 
all the three sources (assuming they are pointlikc) is only 
14.3 ± 4.4, while the number of counts from the entire S- 
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tail is 168 ± 34 (Table 1). Therefore, their contribution 
to the observed X-ray emission from the entire S-tail is 
negligible. 

The PL fits of the outer tails' spectra (see Table 1 
and Fig. 4) indicate that the spectral slope of the S-tail 
does not differ significantly from that of the A-tail (e.g., 
r w 1.5 fits both spectra within the la uncertainties), 
while the N-tail is apparently harder, in contradiction 
to the result found from the XMM-Newton data. Al- 
though the observed number of counts from the N-tail is 
a factor of 1.8 lower than that from the S-tail, their lumi- 
nosities are comparable. Table 1 also suggests that both 
the N-tail and S-tail spectra soften from the bright parts 
toward the extended faint portions, but the statistical 
significance of the softening is low (e.g., Ar w 0.7 ± 0.5 
for the S-tail). The total luminosity of the two tails is 
L w 4.2 x 10 29 d% 50 erg s -1 , in the 0.3-8 keV band, assum- 
ing isotropic emission. This value exceeds the estimate 
derived above from the XMM-Newton data by a factor 
of about 3, but that estimate was obtained for a fraction 
of outer tails, which did not include the bright part of 
the N-tail and included only a small portion of the bright 
part of the S-tail. The average specific intensities in the 
bright parts are about 1.2 x 10 -17 and 2.6 x 10 -17 erg 
cm -2 s -1 arcsec -2 , for the S-tail and N-tail, respectively. 

Similar to the A-tail, the bright and faint portions of 
the S-tail and the faint portion of the N-tail look patchy. 
To quantify the statistical significance of the patchiness, 
we have compared the net counts from brighter and 
fainter regions of equal size for each of these compo- 
nents. For the bright S-tail, faint S-tail, and faint N-tail, 
we found significancies of 2.5 a, 2.2 a, and 2.1 a, respec- 
tively. Therefore, the patchiness is not ruled out, but a 
deeper observation is required to prove it firmly. A sim- 
ilar analysis of the XMM-Newton data does not show a 
statistically significant patchiness, because of large noise 
and poor angular resolution. 

To examine variability of the outer tails, we have com- 
pared the count rates in the area of 462 arcsec 2 of the 
bright part of the N-tail (39.4 ± 8.8 and 5.7 ± 5.0 net 
counts in the observations of 2007 and 2004, respec- 
tively). Accounting for the factor of 1.2 higher sensi- 
tivity of the S3 chip compared to the 13 chip (for the 
spectral parameters derived from the 2007 observation), 
the difference between the count rates is significant at 
the 3.5 a level. On the other hand, the difference of the 
count rates in the 985 arcsec 2 area of the bright part 
of the S-tail, which was detected in both the 2004 and 
2007 Chandra observations (33.0 ± 9.0 and 79.1 ± 13.0 
net counts, respectively), is statistically insignificant. We 
have also looked for variability of the outer tails in the 
XMM-Newton data, but found no statistically significant 
differences between the separate observations because of 
the strong noise. 

As it is natural to assume that the "outer tails" repre- 
sent the sky projection of limb-brightened boundaries of 
a shell, one can expect some X-ray emission from the re- 
gion between the outer tails, in addition to the A-tail. To 
look for this emission, we have inspected two inter-tail re- 
gions of combined area As = 12, 898 arcsec 2 that exclude 
the entire (bright plus faint) A-tail and contain Nt = 886 
counts. Using the background extracted from three 
source-free rectangles around the PWN (Nb = 2222 
counts in the combined area Ab = 29,578 arcsec 2 ), we 



found Ns — —88 ± 36 net counts from the "source". 
Adding the alleged faint portion of the A-tail, which gives 
Nt = 1101 counts in the area As — 15,302 arcsec 2 , we 
found Ns — — 48±41 net counts. Thus, we conclude that 
there is no detectable emission from the region between 
the outer tails, and the "faint portion" of the A-tail is 
likely an illusion (in agreement with our conclusion at the 
end of § 2.2.1). Using the approach outlined by Weisskopf 
et al. (2007), we find the 3 a upper limit of 94 counts in 
the area As = 15, 302 arcsec 2 (99% and 90% upper lim- 
its are 78 and 45 respectively). The corresponding 3 a 
upper limit on the surface brightness, 6.1 x 10 -3 counts 
arcsec -2 , is lower than the average surface brightness of 
the outer tails (e.g., by factors of 4.6 ± 0.9 and 3.4 ± 1.0 
for the "entire" S-tail and N-tail, respectively). 

The XMM-Newton image also does not show de- 
tectable inter-tail emission. Using the same approach, 
we found the 3 a upper limit of 37 counts in the 32" x 64" 
box between the outer tails shown in Figure 1 (99% and 
90% upper limits are 32 and 22 counts, respectively). 
The corresponding 3 a upper limit on surface brightness, 
1.8 x 10 -2 counts arcsec -2 , is lower by factors of 6.2±0.2 
and 2.0 ± 0.4 than the surface brightnesses of the S-tail 
and N-tail, respectively. 

2.2.3. Emission ahead of the pulsar 

An arc-like diffuse emission region, about 5"-7" ahead 
of the pulsar, was reported by P+06 from the 2004 Chan- 
dra data. We have analyzed the data inside a poly- 
gon (area = 86.6 arcsec 2 ) and found 22.8 ± 5.1 source 
counts. Our spectral analysis provides a photon index 
T = 0.9 ± 0.5 and the observed flux F « 1.2 x 10 -14 
erg cm -2 s -1 , in the 0.3-8 keV band, consistent with the 
P+06 estimates. 

Based on these results, the expected number of arc 
counts in the Chandra observation of 2007 7 is 83 ± 24. 
However, although some diffuse emission is seen at that 
site in the 2007 data, there are only 16.3 ± 4.8 source 
counts in the corresponding polygon (at the same dis- 
tance from the pulsar, which has moved 0'.'62 in the sky 
in the 3.5 years), and the shape of the count distribution 
does not resemble an arc (see the upper panel of Fig. 2). 
The spectral slope, T = 1.1 ±0.7, is consistent with that 
of the 2004 arc, but the observed flux, F « 2.5 x 10 -15 erg 
cm -2 s -1 , is a factor of 5 lower. This suggests that the 
emission ahead of the pulsar is variable, but the signifi- 
cance of this variability is not very high (e.g., 2.8 a in the 
difference between the expected and observed counts). 

The alleged arc cannot be seen in the XMM-Newton 
images because it is hidden in the pulsar PSF. However, 
these images show a "streak" ahead of the pulsar, in the 
direction of the proper motion, which is best seen in the 
summed image (bottom panels of Fig. 1). One might 
speculate that this streak is a Geminga PWN element 
(e.g., a forward jet). To check this hypothesis, we ex- 
tracted 760 source counts from the 20" x 50" rectangle 
that includes the streak (shown in the bottom panels of 
Fig. 1) and fit the spectrum with various models. We 
found that the very soft streak spectrum does not fit 
the PL model (x 2 = 4.3 for 44 d.o.f.), but it fits the 
two-component PL+blackbody model (r = 1.8 ± 0.1, 

7 Estimated with the aid of the Chandra PIMMS tool; 
http://asc.harvard.edu/toolkit/pimms.jsp 
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kT = 0.53 ±0.05 keV, xl = 1-1 for 42 d.o.f.) that is con- 
sistent with the pulsar's spectrum (e.g., Kargaltscv et al. 
2005). Therefore, we conclude that the streak is not re- 
lated to the PWN, but it is an artificial spike-like feature 
in the MOS PSF caused by the "spiders" holding the X- 
ray telescopes (see the XMM-Newton Users' Handbook, 
Sec. 3.2.1.1) This conclusion is also supported by a lack 
of any excess above the background at the corresponding 
area in the Chandra images. 

Interestingly, the summed 2004 + 2007 Chandra im- 
age in the bottom panel of Figure 2 shows a hint of a 
short, 10"-15", jetlike structure ahead of the pulsar, with 
its end seemingly connecting to the N-tail. The num- 
ber of counts is, however, too small to conclude whether 
this structure is an accidentally aligned superposition of 
events from the two images or there is indeed a forward 
jet, perhaps bent in the north-northwest direction into 
the N-tail, which we cannot see in the separate images 
because of the scarce statistics. A deeper observation is 
needed to understand the true nature of this and other 
apparent structures in the immediate vicinity of the pul- 
sar. 

3. DISCUSSION 

The Chandra observation of 2007 has confirmed the 
existence of three tail-like features in the Geminga PWN 
image, with comparable luminosities, and allowed us to 
study the PWN in more detail (e.g., to detect the A-tail 
at larger distances from the pulsar and investigate the 
PWN in the immediate vicinity of the pulsar). Moreover, 
it has provided first evidence of variability of the PWN 
elements, in particular, of the A-tail and the emission in 
the pulsar vicinity. 

The observed structure of the Geminga PWN looks 
very unusual. Although the overall appearance of the 
PWN, particularly the alignment of the tails with the 
pulsar's proper motion, leaves no doubts that the PWN 
structure is caused by the supersonic motion of the pul- 
sar in the ISM, none of the other ~20 bowshock-tail 
PWNe detected by Chandra (see KP08) show three dis- 
tinct tails. To interpret the observed structure in terms 
of the PWN models, we should first understand the in- 
trinsic three-dimensional morphology of the X-ray PWN, 
which is by no means obvious. 

At the first glance, the most natural interpretation of 
the PWN elements is that the "outer tails", together 
with the "arc" that apparently connects the tails ahead 
of the pulsar, delineate the limb-brightened boundary 
of the sky projection of an optically thin shell, shaped 
approximately as a paraboloid of revolution, while the 
nearly straight axial tail represents a collimated outflow 
in the direction opposite to that of the pulsar's motion. 
On the other hand, one cannot exclude the possibility 
that the outer tails are, in fact, hose-like structures, such 
as jets confined by their own magnetic fields and bent by 
the head wind of ISM matter. Moreover, one could even 
speculate that the axial "tail" is Dopplcr-boostcd emis- 
sion from a narrow region of a shell formed by material 
flowing with relativistic speeds. Below we will discuss 
these possibilities in more detail. 

3.1. Possible interpretations of the outer tails 

The patchy outer tails of 20"-30" width are seen up to 
~ 3' (<~ 0.2^250 pc) from the pulsar. Their X-ray lumi- 



nosity, L x ~ 4x 10 29 d 2 50 er S s 1 in the 0.3-8 keV band, 
is a fraction of ~ 1.2 x 10 _5 d 2 5 Q of the pulsar's spin-down 
power E, lower than the typical Lx/E <~ 10~ 4 -10~ 3 for 
younger PWNe (KP08). The outer tails' luminosity is 
a factor of ~ 7 lower than the pulsar's magnetosphcric 
luminosity, while the PWN luminosity is usually higher 
than the magnetosphcric luminosity for younger pulsars 
(Kargaltscv et al. 2007). The spectrum of the outer tails 
can be described by the absorbed PL model with T ~ 1, 
which is apparently harder than the typical PWN spec- 
tra. The explanation of the outer tails' properties de- 
pends on the topology of the PW outflow. 

3.1.1. The outer tails are shell boundaries? 

Let us assume that the outer tails are limb-brightened 
shell boundaries and explore the consequences of this 
assumption. First of all, the very fact that the shell 
boundaries are much brighter than the rest of the shell 
image (see § 2.2.2) implies a nonrelativistic speed of the 
bulk outflow along the shell. If it were relativistic, then, 
due to the Dopplcr boosting, the brightest parts of the 
shell image would be not the boundaries but they would 
correspond to the smallest angles between the flow ve- 
locity and the line of sight. This inference restricts the 
number of possible interpretations of the shell. For in- 
stance, the shell cannot be interpreted as synchrotron 
emission from the shocked PW immediately outside the 
(bullet-like) TS surface - not only because the cylindri- 
cal radius of the shell is too large (see P+06), but also 
because the shocked PW is expected to flow with nearly 
relativistic speed, va ow = 0.8c-0.9c, along the lateral TS 
boundary (B+05). As the brightness is proportional to 
[1 — (fflow/c) cos 9]~ r ~ 2 , where 8 is the angle between the 
flow direction and the line of sight, the flow toward the 
observer would be a factor of > 100 brighter than the 
flow in the perpendicular direction, for the photon index 
r w 1. 

One might consider the possibility that the shell emis- 
sion comes from the shocked ISM material heated up to 
X-ray temperatures. In this case, we would associate 
the outer tails with the FS, and the emission mechanism 
with the thermal emission of an optically thin plasma 
rather than synchrotron emission from relativistic elec- 
trons. From the junction conditions at the shock front, 
the expected temperature of the shocked ISM material 
at the head of the bowshock is kT = (3 / '16) u.m p V£ ST — 

0.078yuV2 2 00 keV [where \i is the chemical weight, V200 = 
lp Sr /(200kms _1 ), and adiabatic index 7 a( j = 5/3 is as- 
sumed for the ISM gas], and it should be even lower be- 
hind the pulsar because of the obliqueness of the shock. 
Since the fit of the outer tails spectrum with the mekal 
model gives kT w 5-18 keV (see §2.2.2), the expected 
temperature of the shocked ISM gas looks too low to 
explain the hard spectrum of the tails emission unless 
y psr > 2300M~ 1/2 (£:T/10kcV) 1 / 2 km s" 1 , which would 
imply that the pulsar moves at a small inclination angle 
i with respect to the line of sight, sini < 0.1. However, 
since such a speed is higher than those observed for other 
pulsars, and the PWN appearance can hardly be recon- 
ciled with such small inclination angles, we can discard 
this interpretation. 

It seems more reasonable to assume that the shell is 
formed by the shocked PW flowing immediately inside 
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the CD surface, where the magnetic field is compressed 
(B+05) and the synchrotron radiation is enhanced. How- 
ever, such an interpretation contradicts the available 
PWN models, which predict a nearly relativistic flow 
speed in the outer layers of the synchrotron emitting 
PWN (hence dim boundaries and bright central part) 
because of the Doppler boosting (see above). Indeed, 
the simulated PWN images for i = 90° (see Fig. 4 in 
B+05) show the brightest synchrotron emission from 
the bowshock head region, while throughout the entire 
PWN the brightness decreases from the axis toward the 
CD, in contrast to the observed images. (This is partly 
caused by the assumption that the PWN magnetic field 
is purely toroidal, which reduces the synchrotron inten- 
sity at the PWN boundaries, where the magnetic field is 
parallel to the line of sight. However, even if the mag- 
netic field is completely disordered, the brightness does 
not grow from the axis toward the boundaries; see Fig. 
5 in B+05.) Therefore, this interpretation of the outer 
tails implies that there is a mechanism that decelerates 
the flow. The deceleration can be provided by the shear 
(Kelvin-Hclmholtz) instabilities at the CD, which can 
lead to advection of clumps of the heavier shocked ISM 
material into the shocked PW and slow down the latter 
(e.g., B+05). The diffuse appearance and the patchiness 
of the tails in the high-resolution Chandra images (see 
§ 2.2.2 and Figs. 2 and 3) are consistent with this hy- 
pothesis. It could be verified observationally if the speed 
of the nonuniformities in the outer tails, which should 
form in the process of mass loading, is measured in a 
series of deep observations. 

Another apparent problem with the interpretation of 
the outer tails and the arc as traces of the CD surface 
is the discrepancy between the observed and predicted 
ratios of the CD's cylindrcal radius, tqd, to the distance 
RcB,h ~ 1.3i?TS,h of the CD head from the pulsar. For 
instance, the B+05 model predicts rcD/-RcD,h ~ 3 while 
the observed ratio is > 10 sin i (if we interpret the "arc" 
ahead of the pulsar as the head of the CD surface 8 .) 
This means that either the inclination angle is small, 
i < 15°-20°, or some assumptions of the B+05 model 
are violated. Since the observed PWN shape can hardly 
be reconciled with such small inclination angles (see Fig. 
3 in C+03), we suggest that the discrepancy is caused 
by the assumption that the unshocked PW is isotropic. 
Indeed, if the wind is predominantly equatorial (i.e., con- 
centrated around the plane perpendicular to the pulsar's 
spin axis), and the spin axis is aligned with the direc- 
tion of pulsar motion (as observed in a number of young 
pulsars) , then the lower wind ram pressure ahead of the 
pulsar should result in a smaller distance between the 
pulsar and the TS (and CD) apex. The only models 
of PWNc with anisotropic wind outflow we are aware 
of have been presented by Vigelius et al. (2007). These 
models consider only nonrelativistic flows and, more im- 
portantly, assume zero magnetic field, but they should 
satisfactorily describe the PWN morphology for small 

8 If the PW is anisotropic, the head of the CD surface could 
be much closer to the pulsar than the observed arclike emission. 
In this case, the emission ahead of the pulsar could be a forward 
axial outflow destroyed by the ISM ram pressure. Such an outflow 
should be seen as a diffuse emission of an irregular, variable shape, 
in contrast to the CD surface head that would preserve its arclikc 
shape and could show variations if the ISM is nonuniform. 



values of the PW magnetization parameter (defined as 
the ratio of the Poynting flux to the kinetic energy flux) . 
Although Vigelius et al. did not directly consider the case 
of an equatorial outflow perpendicular to the pulsar ve- 
locity, some of the considered cases (e.g., Fig. 6 in that 
paper) qualitatively confirm our explanation. 

For completeness, we should also mention the interpre- 
tation suggested by C+03, that the outer tails represent 
the synchrotron radiation of the PW in the interstellar 
magnetic field "compressed in the bowshock" by a fac- 
tor of 4 (for an adiabatic shock with 7 ac i = 5/3 and a 
large Mach number), up to ~ 10 /iG. As C+03 assume 
a one-zone shock model, which apparently does not de- 
scribe realistic PWNe, this interpretation may not be 
directly applicable. One may speculate, however, that 
ultrarelativistic electrons from the high-energy tail of the 
electron energy distribution could leak from the shocked 
wind region (incide the CD) into the shocked ISM re- 
gion and generate synchrotron radiation in the interstel- 
lar magnetic field amplified at the FS. We should note, 
however, that only the magnetic field component paral- 
lel to the shock surface is amplified by this mechanism, 
and the amplification becomes insignificant behind the 
pulsar because of the shock obliqueness. Therefore, we 
believe that there is no need to invoke this complicated 
hypothesis as long as the more straightforward explana- 
tion (synchrotron radiation from the region of the CD 
surface) seems viable. 

Finally, we should explain the fact that the outer tails 
are not truly symmetric with respect to the trajectory of 
Geminga on the sky, neither in shape (especially close to 
the pulsar) nor, particularly, in brightness (the S-tail is 
considerably brighter). We could tentatively ascribe this 
asymmetry to nonuniform conditions (density and/or 
temperature) in the ambient medium. The nonunifor- 
mity is supported by the Very Large Array (VLA) and 
Effclsberg radio telescope HI (21 cm line) observations 
of the Geminga field (Giacani et al. 2005). These obser- 
vations have shown the pulsar and its X-ray PWN to be 
in a local minimum of the HI emission, surrounded by 
an open HI shell (an incomplete ring with an average ra- 
dius of 9') that envelopes the southern part of the X-ray 
PWN, with the internal border of the shell close to the 
S-tail. The lack of neutral hydrogen in the vicinity of the 
pulsar can be explained by the ionization caused by the 
pulsar's UV and soft X-ray emission. The openness of 
the shell (no HI emission northwest of the X-ray PWN) 
might imply a higher temperature (and perhaps a lower 
density) of the ISM in that direction. One may specu- 
late that Geminga is crossing a cold ISM cloud and is now 
approaching the cloud's northwest boundary. It remains 
to be understood, however, how the relative brightness 
of the S-tail is connected with the alleged lower tem- 
perature and higher density in that region. We should 
also mention that, based on the HI radio results, one 
could expect an H Q PWN south of the S-tail, associated 
with the FS. However, C+03 report the nondetcction of 
"organized diffuse H Q emission from the X-ray structure 
surrounding Geminga" in a 5 hour observation with the 
VLT-ANTU telescope. Their Figure 2 shows an appar- 
ent filament at the outer border of the S-tail, but it is not 
immediately clear whether or not this feature is related 
to the PWN. 

To conclude, if the outer tails represent the limb- 
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brightened boundaries of the sky projection of a shell, 
this shell is most likely the synchrotron radiation from 
the region of interaction of the shocked PW and shocked 
ISM material, where the wind flow is decelerated to non- 
relativistic velocities by the shear instability, which im- 
plies mass loading. The shape of the shell is somewhat 
different from the shape of the CD surface in the available 
numerical PWN models, perhaps because the models do 
not include the mass loading and proper anisotropy of 
the unshocked PW. 

3.1.2. The outer tails are bent polar outflows? 

The hypothesis that the outer "tails" represent the 
boundary of a shell is not the only possible explanation. 
In particular, as no emission is seen between the outer 
tails except for the A-tail (see §2.2.2), we cannot exclude 
the possibility that the outer tails are in fact some colli- 
mated flows emanating from the pulsar magnetosphere, 
such as two jets aligned with the pulsar's spin axis near 
the pulsar and bent by the ram pressure at larger dis- 
tances. This interpretation implies a large angle 9 be- 
tween the spin axis and pulsar's velocity (9 > 60°-70°, 
as follows from Figs. 2 and 3) and a sufficiently large 
angle £ between the spin axis and the line of sight 9 . 

Chandra observations have shown that jets emanating 
along the spin axes are ubiquitous among young PWNe 
(see, e.g., Weisskopf et al. 2000; Pavlov et al. 2003), and 
the spin axis is often approximately aligned with the pul- 
sar velocity direction (Ng & Romani 2007). The mecha- 
nisms of jet formation and collimation are currently not 
certain. In the scenario discussed by Benford (1984), a 
fraction of electrons created in the vacuum gaps above 
the magnetic poles and accelerated along the open mag- 
netic field lines is deflected toward the spin axis and 
forms a beam collimated by its own toroidal magnetic 
field. Another mechanism of axial outflow formation has 
been discussed by Komissarov & Lyubarsky (2004), who 
assume that the outflow is originally equatorial and show 
that the magnetic hoop stress can stop the outflow in the 
surface layers of the equatorial disk and redirect it into 
magnetically confined polar jets. 

If the pulsar were not moving with respect to the am- 
bient medium, the jet matter would keep flowing along 
the spin axis until the jet is destroyed by the interaction 
with the medium. The ram pressure exerted onto the 
jets of a moving pulsar can bend the jets in the direction 
opposite to that of the pulsar's motion, so that the jets 
are seen as two tails behind the moving pulsar 10 . 

This scenario allows one to explain the observed asym- 
metry of the Geminga's outer tails (see §2.2.2). The 
asymmetry can be associated with the large (but differ- 
ent from 90°) angle 9 between the spin axis and the pul- 
sar's velocity. At such an orientation the angles between 
the ram pressure direction and the jet matter velocity 
directions are acute and obtuse for the southeastern and 
northwestern jets, respectively, which means that bend- 
ing the northwestern jet is more difficult. An additional 
reason for the asymmetry might be a deviation of the 
pulsar's magnetic field geometry from an ideal centered 

9 The large value of £ supports the outer gap interpretation of 
Geminga's 7-ray emission (Romani & Watters 2010). 

10 A similar model has been discussed by Heinz et al. (2008) for 
microquasars. 



dipole, which would lead to different structures of the 
magnetic field at the two poles and different properties 
of the two jets. Different brightness of the jets, espe- 
cially in the pulsar vicinity, might be caused by Doppler 
boosting (if the angle ( between the spin axis and line of 
sight is different from 90°), but it is hard to estimate the 
Doppler factor and to infer the angles with the current 
noisy data. 

As the bent polar outflow interpretation of the outer 
tails requires a large value of the angle 9, while 9 is ap- 
parently small for most pulsars, this interpretation im- 
plies that the outer tails are a rare phenomenon, in agree- 
ment with PWN observations that have not shown such 
tails in other PWNe. It, however, remains to be under- 
stood why 9 is so different for Geminga. A theoretical 
study of the expected distribution of this angle using the 
physics of the neutron star birth is required to confirm 
this explanation. 

3.2. The nature of the axial tail 

The straight, patchy A-tail is seen up to 50" (0.06c?250 
pc) from the pulsar. Its surface brightness is not only 
nonuniform but also variable, as we see from the com- 
parison of the 2007 and 2004 data. Assuming a nearly 
isotropic emission (which, rigorously speaking, implies 
a nonrelativistic flow), the A-tail luminosity, Lx ~ (1- 
2) x 10 29 ^ 50 erg s - * in the 0.3-8 keV band, is - (3- 
6) x 1 0~ 6 <i25o °f Geminga's spin-down power E. The A- 
tail luminosity is ~ 0.03-0.06 of the nonthermal (magne- 
tospheric) luminosity of the Geminga pulsar in the same 
energy range and is a factor of 2-4 lower than the total 
luminosity of the outer tails. 

There are three conceivable explanations of the A-tail: 
a jet emanating from the pulsar magnetosphere in the 
direction opposite to the pulsar velocity, a tail part of 
the bowshock-tail PWN created by the supesonic mo- 
tion of the pulsar, and a Doppler-boosted image of a shell 
into which a fraction of the relativistic PW is directed. 
We will discuss these interpretations below, taking into 
account their connection with the above-discussed inter- 
pretations of the outer tails. 

3.2.1. The axial tail is a pulsar jet? 

The interpretation of the A-tail as a pulsar jet, sug- 
gested by P+06, is consistent with only one of the above- 
discussed interpretations of the outer tails, namely, the 
hypothesis that the outer tails represent a boundary of 
a shell (e.g., an equatorial outflow bent by the ram pres- 
sure). As pulsar jets emerge along the pulsar's spin axis, 
and the A-tail is aligned with the pulsar's trajectory in 
the sky, the jet interpretation of the A-tail implies that 
the spin axis is likely aligned with the pulsar's velocity. 
This suggests that the "natal kick" of the Geminga pul- 
sar was directed along the spin axis, which has important 
implications for the mechanisms of supernova explosion 
and neutron star formation (e.g., Ng & Romani 2007). 

The lack of a clear (counter )jet ahead of the pulsar 
(see, however, the note at the end of § 2.2.3) could be 
explained by Doppler boosting (the approaching jet is 
brighter than the receding counterjet, assuming the jet 
material flows with nearly relativistic velocities). Alter- 
natively, the counterjet can be partially or fully destroyed 
by the ISM ram pressure, or the outflows in the opposite 
directions may be intrinsically different. 
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As described in §2.2.1, the surface brightness is 
distributed nonuniformly along the A-tail, with some 
"blobs" seen at different positions in the images of 2004 
and 2007, and a "gap" between the pulsar and the begin- 
ning of the A-tail. This means that there are regions of 
the enhanced magnetic field and/or higher density along 
the A-tail, which might be caused by discrete ejections 
from the magnetosphere, or they could be manifestations 
of some instabilities (e.g., the sausage instability, as dis- 
cussed by Pavlov et al. 2003 for the outer jet of the Vela 
PWN) or internal shocks in the jet flow. Particularly in- 
teresting is the brightest blob C seen at the apparent end 
of the A-tail in the image of 2007 (see Fig. 2). Although a 
background K star could contribute to the blob C emis- 
sion (see §2.2.1), our analysis suggests that the star's 
contribution is not dominant and blob C could be asso- 
ciated with the jet's termination shock. To understand 
the nature of the blobs, it would be important to study 
their evolution in a series of deep observations, which, in 
particular, would help estimate the flow speed in the jet. 

Using the observed diameter of the alleged jet and the 
estimate for the energy injection rate, P+06 estimate 
the jet's magnetic field: B }et ~ 100 fiG. For such a mag- 
netic field the expected jet length is Zj ct ~ fjctT S yn ~ 
0.6 (ujet/0.5c) (Bjet/lOO/uG)- 3 / 2 pc, where Wjot is the 
bulk flow velocity in the jet, and T syn is the synchrotron 
cooling time. The jet length estimated from the initial 
bright portion of the axial tail, (j ct ~ 0.06 c^so/ sini pc, is 
much smaller than this value unless Vj et is much smaller 
than 0.5c and/or sinz is small, which seems unlikely To 
explain this contradiction, P+06 speculate that the jet 
becomes uncollimatcd or destroyed well before it radiates 
its entire internal energy. To check such speculations and 
test the pulsar jet interpretation, deeper observations are 
required. 

3.2.2. The axial tail is a shocked pulsar wind? 

A tail-like structure similar to the observed A-tail 
could form behind the supersonically moving pulsar due 
to the collimation of the shocked PW by the ram pres- 
sure (B+05; Romanova et al. 2005). For instance, if we 
assume that the outer tails and the possible arc ahead 
of the pulsar delineate the CD surface (see §3.1.1), then 
the A-tail might be the shocked PW immediately out- 
side (and perhaps behind) the bullet-like TS. P+06 have 
shown that this interpretation is not quantitatively con- 
sistent with the available simulations of bowshock-tail 
PWNc, but those simulations do not take into account 
the intrinsic anisotropy of the PW. In addition, it would 
be difficult to explain the presence of the blobs and the 
variability of the A-tail in the framework of this inter- 
pretation. Therefore, we consider this interpretation un- 
likely. 

If the outer tails are bent polar outflows (§ 3.1.2), then 
the ram-pressure confined shocked PW would be the 
only possible explanation for the A-tail. A more de- 
tailed interpretation of the A-tail would depend on the 
PW model. For instance, if the PW were intrinsically 
isotropic, then the A-tail might be interpreted as origi- 
nating from the shocked PW "sheath" immediately out- 
side the TS, and the observed width of the tail would 
imply the distances i?TS.h ~ 3" and i?cD.h ~ 4" of the 
TS and CD heads from the pulsar. Therefore, we would 
expect a bright arc (brighter than any part of the ax- 



ial tail) « (3"-4")sini ahead of the pulsar. No such a 
bright arc is seen in the images, but this does not nec- 
essarily rule out the TS origin of the A-tail because the 
TS head could be closer to the pulsar and hidden within 
the pulsar image if the PW is anisotropic (P+06). In 
this interpretation, however, we have to assume that no 
emission is seen from the CD surface region, which looks 
somewhat unnatural. 

In the bent polar outflow interpretation of the outer 
tails, one could also assume that the A-tail tail is associ- 
ated with the CD-confined cylindrical region behind the 
(unresolved) TS. In this case, for an isotropic PW, we 
would expect a bright arc ~ l"sini ahead of the pulsar, 
which can easily be hidden within the pulsar image. In 
the framework of this interpretation, the non-uniformity 
and variability of the A-tail could be explained by shear 
instabilities at the CD surface (cf. §3.1.1), which could 
also decelerate the flow, so that the blobs' velocity would 
be lower than that in the jet interpretation of the A-tail. 
We should note, however, that in such interpretations the 
polar outflows are more luminous than the ram-pressure 
confined tail (perhaps an equatorial outflow), which has 
not been observed for any other pulsar. 

3.2.3. The axial tail is a Doppler-boosted shell 
projection ? 

As we have mentioned above, the image of a shell 
formed by material outflowing with relativistic speeds 
may be strongly affected by Doppler boosting, which 
brightens those parts of the shell where the angle 9 be- 
tween the bulk flow velocity and the line of sight is the 
smallest. For instance, a conical shell, in which the ma- 
terial flows from the cone vertex at the pulsar position, 
would look like a straight strip 11 corresponding to the 
minimum 9. If the material flows from the head of a 
paraboloid-like shell (e.g., the CD surface), the observer 
would see a shorter strip, detached from the paraboloid 
head in general case. In principle, one could imagine 
that the A-tail is such a projection of the shell formed 
by the shocked PW that flows out with relativistic ve- 
locities between the TS and CD, while the outer tails 
are bent polar outflows (jets) . In this interpretation, the 
true transverse radius of the shell would be larger than 
the observed width of the A-tail. The blobs in the A- 
tail might be some local instabilities in the relativistic 
flow, which would likely move with relativistic bulk-flow 
velocities. Therefore, it would be important to examine 
the blob motion in future observations. 

4. SUMMARY 

The new Chandra and XMM-Newton observations of 
the Geminga PWN have confirmed that it has three 
tail-like components, unlike any other detected PWN. 
The new observations have allowed us to image the 
tails at larger distances from the pulsar and establish 
their patchy structure. Comparing the new and previ- 
ous Chandra observations, we have found indications of 
PWN variability, especially in the axial tail and the emis- 
sion ahead of the pulsar. In particular, we found up to 

11 A similar model has been discussed by Radhakrishnan & 
Deshpande (2001) who interpret the "inner jets" of the Vela PWN 
as a projection of the rotating beams of relativistic particles emit- 
ted along the magnetic axes. 
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three blobs in the axial tail, at different positions in 2004 
and 2007. 

Similar to other X-ray PWNe, the Geminga PWN is 
due to synchrotron radiation of shocked PW comprised 
of relativistic particles. Based on the new and old obser- 
vations, we have proposed several competing interpreta- 
tions of the PWN structure. Very likely, the outer tails 
delineate a limb-brightened boundary of a shell-like re- 
gion of interaction of the shocked PW and shocked ISM, 
while the axial tail is a pulsar jet along the spin axis 
aligned with the pulsar's trajectory. Such an interpreta- 
tion implies a nonrelativistic speed of the bulk outflow 
along the shell, possibly decelerated by the shear insta- 
bility and mass loading. Alternatively, the outer tailis 
could be polar outflows from the pulsar magnetosphere 
(e.g., pulsar jets along the spin axis), bent by the ISM 
ram pressure, in which case the axial tail could be a 
shocked PW (e.g., an equatorial outflow) collimatcd by 
the ISM ram pressure exerted on the supersonically mov- 
ing PWN. 

To discriminate between various interpretation of the 
observed PWN, a series of carefully designed Chandra 
observations is required. In particular, such observations 
should allow one to measure the speeds of the bulk flows 
in the tails, which would distinguish fast jets from ram- 
pressure-confined pulsar winds slowed down by the in- 



teraction with the ambient ISM. Also, such observations 
should be deep enough to establish the true morphology 
of the emission in the immediate vicinity of the pulsar. 
For instance, if a deeper observation convincingly shows 
that there is an arc ahead of the pulsar connecting the 
two outer tails, then the bending axial outflows scenario 
will be ruled out. If, however, we see two straight tails 
originating from the pulsar in a direction inclined to the 
pulsar's velocity direction, then the tails can be inter- 
preted as bent jets. In addition, the detailed modeling of 
anisotropic magnetic PW from a high-speed pulsar will 
also be extremely useful to properly interpret the obser- 
vational data. 
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New X-ray Observations of the Geminga PWN 



TABLE 1 

Properties of the PWN elements from the Chandra data of 2007 and 2004. 
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Note. — The initial bright and extended faint portions of the three PWN tails in the data of 2007 arc 
shown in Fig. 3. The A-tail, N-tail, and S-tail stand for the axial tail, the northern tail, and the southern 
tail, respectively. The spectra of the PWN elements were fitted with the absorbed power-law model (photon 
index T) at a fixed hydrogen column density TVh = 1-1 X 10 20 cm -2 . The errors are given at the 1<t level, 
for one interesting parameter. 

a Area of the PWN element in arcsec 2 . 
b Background-subtractcd counts in PWN elements. 

c Observed (absorbed) flux in the 0.3—8 keV band, in units of 10 — 14 erg cm -2 s — 1 , estimated with XSPEC. 
d Unabsorbcd luminosity in the 0.3-8 keV band, assuming isotropic emission and distance of 250 pc, in units 
of 10 29 erg s -1 . 

e Blob C at the end of the brighter portion of the axial tail is included. 
f Blob C (3" radius) is excluded from the data. 
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Fig. 1.— XMM-Newton MOS1+MOS2 images (5' X 5') of the Geminga pulsar and its PWN in the 0.5-8 keV band. The top, middle 
and bottom panels correspond to the observation of 2004-2007 (64 ks total scientific exposure), 2002 (78 ks), and 2002-2007 (142 ks), 
respectively. The images in the left panels are binned in 3" X 3" pixels, while the images in the right panels are additionally smoothed with 
a 18" FWHM Gaussian. The ellipses (120" X 40") show the regions for which the spectra and fluxes were measured, while the 64" X 32" 
rectangle between the ellipses was used for estimating the upper limit on the surface brightness between the outer tails (see §2.2.2). The 
50" X 20" rectangle ahead of the pulsar was used to measure the spectrum of the "streak" (see § 2.2.3). The 5" radius circle in the bottom 
left panel is centered at the position of blob C found in the Chandra observation of 2007 (see Fig. 2 and § 2.2.1). The source north-northwest 
of the pulsar is an 11-th magnitude K star (C+03). 
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Fig. 2. — lf6 X If 6 images of the Geminga pulsar and its PWN from the Chandra observations of 2007 and 2004 (top and middle panels, 
respectively), and the combined image (bottom), in the 0.3-8 keV range. The original data were binned in 1" X 1" pixels and smoothed 
with a 2" FWHM Gaussian. The arrows show the pulsar's proper motion, whereas the dashed lines in the middle and bottom panels 
indicate the sky region for which the exposure was reduced by the telescope dithering (from 18.8 to ks, in the 2004 observation). Three 
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Fig. 4. — Confidence contours for the PL normalization versus photon index (at the 68.3% and 95.4% levels, for two interesting parameters) 
for bright portions of the PWN tails from the Chandra 2007 data. The hydrogen column density is fixed at Nh = 1.1 x 10 20 cm -2 . The 
dotted curves correspond to constant unabsorbed flux values (depicted near the curves, in units of 10 — 14 erg cm~ 2 s _1 ). 
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Fig. 5. — Count rate profiles along the axial tail in the Chandra data. The net count rates are calculated in square boxes of equal size 
(59.5 arcsec 2 ), shown in the inset. The dashed (black) and solid (red) points (with la errors) are for the 2004 and 2007 data, respectively. 
The count rates and their errors for the 2007 data are multiplied by a factor of 1.6 to account for the different sensitivities of the 13 and 
S3 chips. 



16 



Pavlov, Bhattacharyya, Zavlin 




6 h 34 m 03?Q 33 m 56?0 49?0 42?0 

R.A. 

Fig. 6. — 5.' 4 X 5'4 image of the Geminga pulsar and its PWN from the Chandra data of 2007 binned in 1" X 1" pixels and smoothed 
with an 8" FWHM Gaussian, with overlaid brightness contours from the XMM-Newton image shown in Fig. 1. The green crosses with 
numbers indicate the positions of optical stars projected onto the PWN elements. 



